Isospin dependence of projectile-like fragment production at intermediate energies 
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The cross sections of fragments produced in 140 A MeV 40 ' 48 Ca + 9 Be and 58 ' 64 Ni + 9 Be reactions 
are calculated by the statistical abration-ablation(SAA) model and compared to the experimental 
results measured at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State 
University. The fragment isotopic and isotonic cross section distributions of 40 Ca and 48 Ca, 58 Ni and 
64 Ni, 40 Ca and 58 Ni, and 48 Ca and 64 Ni are compared and the isospin dependence of the projectile 
fragmentation is studied. It is found that the isospin dependence decreases and disappears in the 
central collisions. The shapes of the fragment isotopic and isotonic cross section distributions are 
found to be very similar for symmetric projectile nuclei. The shapes of the fragment isotopic and 
isotonic distributions of different asymmetric projectiles produced in peripheral reactions are found 
very similar. The similarity of the distributions are related to the similar proton and neutron density 
distributions inside the nucleus in framework of the SAA model. 

PACS numbers: 25.70.Mn, 21.65.Cd 
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I. INTRODUCTION 

Projectile fragmentation is a well-established technique 
for the production of rare isotope beams used by many 
radioactive ion-beam facilities around the world. The 
process of projectile fragmentation has been studied ex- 
tensively to investigate the reaction mechanisms in heavy 
ion collisions at intermediate and high energies. Under- 
standing the physics of projectile fragmentation is impor- 
tant not only for rare-isotope beam production purposes 
but also for the fundamental nuclear physics processes 
involved in nuclear collisions [l|, |H, H, 0, H, H, 0] • 

Isospin effect is the phenomenon induced by the isospin 
degree of freedom in heavy- ion collisions. Isospin effects 
of various physical phenomena, such as multifragmen- 
tation, collective flow, preequilibrium nucleon emission, 
etc., have been reported [J i, ©, M, El 03, 03, G|. 
These studies have shown that isospin effect exists in nu- 
clear reactions induced by exotic nuclei but it may dis- 
appear under some conditions. In projectile fragmen- 
tation reactions, the yields of neutron-rich nuclei from 
fragmentation of neutron-rich projectiles will be larger 
than those from stable nuclei; this is one isospin effect in 
fragment production[T^. [TEl fill [L7| . But the difference 
becomes smaller with the increase of the charge difference 
between the fragment and the projectile. And it may dis- 
appear at last. In Refs. [l5| and [l6j], the fragment iso- 
topic distributions of 60 A MeV 16 ' 18 + 9 Be, 36 > 40 Ar + 
9 Be and 40 ' 48 Ca + 9 Be were presented by Fang et al. It 
was found that the peak position of isotopic distributions 
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from stable nucleus induced reactions has a shift toward 
the neutron-rich side as compared to that from neutron- 
rich nucleus induced reactions. The shift becomes larger 
when the difference between neutron numbers of the two 
projectiles is bigger. This phenomena is called the isospin 
effect in the projectile fragmentation. But the isospin ef- 
fect of fragmentation reaction on the isotopic distribution 
decreases with the increase of the atomic number differ- 
ence and disappears at last. 

In the campaign of four projectile fragmentation ex- 
periments carried out at the National Superconduct- 
ing Cyclotron Laboratory (NSCL) at Michigan State 
University [3, eight different reaction systems and more 
than 1400 fragment cross sections have been measured. 
Reactions of primary 140 A MeV of 40 Ca, 48 Ca, 58 Ni, 
and 64 Ni colliding with light target 9 Be and heavy tar- 
get 181 Ta were measured. Extensive study of the projec- 
tile fragmentation reactions using the EPAX code|18|. 
the macroscopic-microscopic heavy ion phase space ex- 
ploration (HIPSE) model, and the fully microscopic anti- 
symmetrized molecular dynamics (AMD) model has been 
carried out [191]. With no variation of the model parame- 
ters, a reasonable agreement between the predictions and 
experimental data has been reached [l9j. Comprehensive 
fragment cross sections of 40 Ca_, 48 Ca, 58 Ni and 64 Ni were 
presented in Refs. [l|| and [lj|. From the detailed cross 
section of fragments presented, it is possible to investi- 
gate the isospin effect of projectile fragmentation induced 
by the symmetric and asymmetric (neutron-rich) nuclei. 

One of the methods used to calculate the cross sec- 
tion of fragments produced in projectile fragmentation 
is the statistical abrasion ablation (SAA) model. Com- 
pared to the HIPSE model and the AMD model, the cal- 
culation of the SAA model is simple and can reproduce 
the experimental results of heavy ion collision at inter- 
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FIG. 1: (Color online) The cross sections of fragments produced in the 140 A MeV 40,48 Ca + 9 Be reactions. The open circles 
denote the fragment isotopic a distributions, and the filled ones denote the fragment isotonic a distributions of the 40 Ca + 9 Be 
reaction. The x axis represents the numbers of protons removed(AZ = Zp ro j — Zf ra g) from the projectile for the isotones, and 
the numbers of neutrons removed(AN = Np ro j — Nf ra g) from the projectile for the isotopes in the pad. The open triangles 

denote the fragment isotopic a distributions and the solid triangles denote the fragment isotonic a distributions of 48 Ca + 9 Be 
reaction. The lines are the results of the SAA model calculations. 



mediate energy well[H d M EH El- In this article ' 
the cross section of fragments produced in 140 A MeV 
40 ' 48 Ca + 9 Be and 58 > 64 Ni + 9 Be reactions are calculated 
within the framework of the SAA model and compared to 
the NSCL experimental data [19]. The fragment isotopic 
and isotonic cross section distributions between the four 
reactions is compared to study the isospin dependence 
of projectile-like fragmentation for symmetric and asym- 
metric nuclei. 



"participants" and "spectators." The participants in an 
overlapping region between the projectile and the target 
interact strongly while the spectators are left to move al- 
most without being disturbed 25]. In the second evapora- 
tion stage, the system reorganizes because of excitation, 
which means that it is deexcited and thermalized by the 
cascade evaporation of light particles. After the deexcita- 
tion, the results of the final fragment, which are compa- 
rable to the experimental data, can be obtained. The de- 
tails of the SAA model can be found in Refs. [HE3,E1- 



II. THE SAA MODEL 

The SAA model was developed by Brohm and 
Schmidt to describe the peripheral nuclear collisions at 
high energies in a picture of quasi-free nucleon-nuclcon 
collisions [23[. It was modified by Fang and Zhong 
et al. to study the heavy ion collisions at intermedi- 
ate energy [H |l3, ES EJ, E3, El- The reactions in the 
SAA model are described as a two-step process. The ini- 
tial stage can be described by a Glauber-type model as 



In the SAA model, the colliding nuclei are described to 
be composed of parallel tubes orienting along the beam 
direction. Neglecting the transverse motion, the collision 
is described by independent interactions of tube pairs. 
Assuming a binomial distribution for the absorbed pro- 
jectile neutrons and protons in the interaction of a spe- 
cific pair of tubes, the distributions of the total abraded 
neutrons and protons can be determined. For an in- 
finitesimal tube in the projectile, the transmission prob- 
abilities for neutrons (protons) at a given impact param- 
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FIG. 2: (Color online) The cross sections of fragments produced in the 140 A MeV 58,64 Ni + 9 Be reactions. The open circles 
denote the fragment isotopic a distributions, and the solid ones denote the fragment isotonic a distributions of the 58 Ni+ 9 Be 
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from the projectile for the isotones, 



and the numbers of neutrons removed(AN=Np ro j— Nf ra g) from the projectile for the isotopes in the pad. The open triangles 



denote the fragment isotopic a distributions and the solid triangles denote the fragment isotonic a distributions of the 
9 Be reaction. The lines denote the results of the SAA model calculations. 
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eter b are calculated by 

tfc(s-b) =exp{-[^(s-bK fe + -D^(s-b)(T pfe ]} (1) 

where D T is the nuclear-density distribution of the tar- 
get intergrated along the beam direction and normalized 
by jd 2 sDl = N T and J d 2 sD^ = Z T . N T and Z T 
refer to the neutron and proton numbers of the target, 
respectively. The vectors s and b are denned in the plane 
perpendicular to the beam, a^k is the free space nucleon- 
nucleon cross section (k', k = n for neutron and k',k = p 
for proton). The average absorbed mass in the limit to 
infinitesimal tubes at a given b is 



< AA(b) >- 



d 2 sDl(s)[l-t n (s-b)} 



dW(s)[l 



t p (s - b)] 



(2) 



The excitation energy of projectile spectator is estimated 
by the simple relation of E* = 13.3 < AA(b) > MeV, 
where 13.3 is the mean excitation ener gy due to an 
abraded nucleon from the initial projectile[26|. 



The production cross section for a specific isotope can 
be calculated from 



a(AN,AZ)= J d 2 bP(AN,b)P(AZ,b), 



(3) 



where P(AN, b) and P(AZ, b) are the probability distri- 
butions for the abraded neutrons and protons at a given 
impact parameter 6, respectively. 

The a of fragments produced in 140 A MeV 40 > 48 Ca 
+ 9 Be and 58 ' 64 Ni + 9 Be reactions are calculated within 
the framework of the SAA model. In the calculations, the 
free space nucleon- nucleon cross sections are adopted [27|. 
The proton and neutron density distributions are as- 
sumed to be the Fermi type, 



Pi(r) 



Pi 



l + ex P (^)' 



n,p 



(4) 



where p° is the normalization constant that ensures that 
the integration of the density distribution equals the 
number of neutrons (i = n) or protons [i = p); ti is the 
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diffuseness parameter, and C, is half the density radius 
of the neutron or proton density distribution. 

III. RESULTS AND DISCUSSION 

In Fig. [U the a of fragments produced in 140 A MeV 
40 ' 48 Ca + 9 Be are plotted. In Fig. [2J the a of fragments 
produced in 58,64 Ni+ 9 Be reactions are plotted. The x 
axis is not the usually used mass number but the number 
of removed neutrons(AiV = -^p ro j — -^frag) ^ or isotopes 
or the number of removed protons(AZ = -Zp ro j — ^frag) 
for isotones from the projectile, which not only can re- 
flect the fine information of the reactions, but also makes 
it convenient to compare between the fragment isotopic 
and isotonic distributions of different projectiles. From 
Figs. Q] and [2j it can be seen that the SAA model re- 
produces the experimental data quite well not only for 
the stable nucleus projectile but also for the neutron-rich 
nucleus projectile. 

In Refs. and the fragment isotopic distribu- 
tions of different projectiles are compared. The isospin 
effect and its disappearance in the projectile fragmenta- 
tion is found. There is a shift from the fragment isotopic 
a distributions (fragment isotopic distribution) of 40 Ca 
to that of 48 Ca. There is also a shift from the isotonic 
a distributions of 40 Ca to that of 48 Ca. It is the same 
phenomena as observed in Refs. [l5[ and [IB]; i.e., for the 
neutron-rich projectile, neutrons are more easily removed 
than for the symmetric nucleus projectile. 

For 40 Ca, in Fig. [TJ the shapes of the fragment isotopic 
and isotonic distributions are similar but there is a shift 
from the isotopic to the isotonic distributions. The shift 
indicates the isospin dependence of the fragment produc- 
tion. The shift becomes smaller when AZ of fragment 
isotopes and AN of fragment isotones increases, which 
means that central collisions begin to dominate the re- 
actions. In other words, the isospin effect decreases in 
central collisions. The fragment isotopic distributions of 
48 Ca are wider than those of 40 Ca. It is easy to under- 
stand because 48 Ca is a neutron-rich nucleus and has a 
neutron skin structure, from which neutrons can be re- 
moved more easily. The shift between the fragment iso- 
topic and isotonic distributions of 48 Ca is in the reverse 
direction of that of 40 Ca and become wider compared to 
those of 40 Ca, which indicates that the isospin effect in 
48 Ca projectile fragmentation does not decrease but is 
enhanced in central collisions. 

The shapes of fragment isotopic and isotonic distribu- 
tions of 58 Ni are similar to those of 40 Ca. The fragment 
isotopic and isotonic distributions of 64 Ni are similar to 
those of 48 Ca. What interested us most in Fig. [2] is that 
the fragment isotopic and isotonic distributions overlap 
when AZ of isotopes and A7V of isotones are more than 
9. It is the direct evidence of the disappearance of the 
isospin dependence of the projectile fragmentation. Com- 
pared to that of 40 Ca, it also means there is less isospin 
dependence in the projectile fragmentation of 58 Ni. 



From the investigations of the fragment isotopic and 
isotonic distributions of 40 Ca and 58 Ni, evidence for the 
decrease and disappearance of the isospin effect in projec- 
tile fragmentation is found; i.e., the shift from the frag- 
ment isotopic distributions to the fragment isotonic dis- 
tributions of 40 Ca and 58 Ni becomes narrower and the 
fragment isotopic and isotonic distributions of 58 Ni over- 
lap. 

There is little difference between the proton and neu- 
tron density distributions inside the symmetric nucleus. 
For 40 Ca and 58 Ni, with N/Z equal to 1.0 and 1.071, 
respectively, the proton and neutron density distribu- 
tions should be very similar except the size of 58 Ni is 
larger than that of 40 Ca. There are some reasons to ex- 
pect similar fragment isotopic and isotonic distributions 
of their projectile fragmentation. In Fig. [31 the a val- 
ues of fragments produced in 40 Ca/ 58 Ni + 9 Be reactions 
are plotted. Though the fragment isotopic distribution of 
40 Ca is narrower than that of 58 Ni, the difference between 
the fragment isotopic distributions of 40 Ca and 58 Ni is 
small. The same behavior is exhibited in the fragment 
isotonic distributions. Because the difference between 
40 Ca and 58 Ni mass numbers is not very big, it can be 
concluded that, for symmetric projectile nuclei with sim- 
ilar mass numbers, the fragment isotopic and isotonic 
distributions should be similar. Studying these distribu- 
tions more carefully, it can be seen that the a values of 
fragments produced in peripheral reactions of Ca and 
58 Ni have less difference than those in central collisions. 
For the AZ > 3 isotopes, the a values of fragments of 
58 Ni are higher than those of 40 Ca on the right sides of 
the distributions while the left sides of the distributions 
show very little difference. For the AN > 3 isotones, the 
a values of fragments of 58 Ni are higher than those of 
40 Ca on the left sides of the distributions while the right 
sides of the distributions show very little difference. 

In Fig. [H the u values of fragments produced in 48 Ca 
/ 64 Ni + 9 Be reactions are plotted. There is very little 
difference between the fragment isotopic distributions of 
48 Ca and 64 Ni in peripheral reactions. For the AZ < 3 
isotopes, though the right sides (more nucleons removed) 
of the isotopic distributions have a shift from the frag- 
ments of 48 Ca to 64 Ni, the left sides (little nucleons re- 
moved) of the isotopic distributions show very little dif- 
ference. The AN < 3 isotones of 48 Ca and 64 Ni overlap. 
The fragments on the left side of the AZ < 3 isotopes 
and the AN < 3 isotones are the productions of most pe- 
ripheral reactions. In Fig. [J] the shapes of the AZ > 3 
isotopic distributions of 48 Ca and 64 Ni are very similar 
but the isotopic distributions of 48 Ca shift to those of 
64 Ni and the right sides of the distributions overlap grad- 
ually. Compared to the left side of the distribution, the 
right side of the distribution means that the central colli- 
sions happen most. For 48 Ca and 64 Ni, the distributions 
of AZ < 5 fragment isotones, which are the productions 
in peripheral reactions, show very little difference. 

By investigating the fragment distributions of the four 
projectiles, the isospin effect in the projectile fragmen- 
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FIG. 3: (Color online) The cross sections of fragments produced in the 140 A MeV 40 Ca + 9 Be and 58 Ni + 9 Be reactions. 
The open circles denote the fragment isotopic a distributions, and the solid ones denote the fragment isotonic a distributions 
of the 40 Ca + 9 Be reaction. The open triangles denote the fragment isotopic a distributions and the solid triangles denote the 
fragment a distributions of the 58 Ni+ 9 Be reaction. The lines denote the results of the SAA model calculations. 



tation is exhibited. By comparing the fragment isotopic 
and isotonic distributions of 40 Ca and 58 Ni, it has been 
found that the isospin dependence of the projectile frag- 
mentation decreases and disappears. The similarities of 
the fragment isotopic and isotonic distributions from pro- 
jectile fragmentation are discovered not only in symmet- 
ric nuclei but also in asymmetric nuclei. It is meaningful 
to discover the similarities of the fragment isotopic and 
isotonic distributions because it help us to estimate the 
production of fragments in projectile fragmentation at 
intermediate energy. But why do these phenomena hap- 
pen? 

According to the SAA model, the number of protons 
or neutrons removed in a specific pipe in collision is de- 
termined by the nucleon-nucleon cross section and the 
densities of the protons and neutrons. Because the pro- 
jectile fragmentation of 40,48 Ca and 58 ' 64 Ni occur under 
the same experimental conditions, the deexcitation in the 
reactions is the same and there is no need to consider the 
target effect in the reactions. Some qualitative explana- 
tions of the phenomena observed above can be obtained 
from the SAA model. 

For symmetric nuclei, as have been referred, their pro- 



ton and neutron density distributions are similar. Ac- 
cording to the SAA model, the fragment isotopic and iso- 
tonic distributions of projectile fragmentation should be 
similar and the difference should reflect the difference of 
the proton and neutron distributions. The shift from the 
fragment isotopic distribution to the isotonic distribution 
in peripheral reactions of 40 Ca reveals the difference of 
the proton and neutron density distributions in the sur- 
face region. It should be kept in mind that even in the 
central collisions the surface effect of nucleus also exists 
because in a specific pipe involved in collisions, the den- 
sity distribution is the mixed effect of the surface and the 
core. But the surface effect decreases in the central colli- 
sion since the density of core will wash out some density 
difference in the surface region. It make the shift from 
the fragment isotopic to isotonic distributions become 
smaller. For 58 Ni, the surface effect become unimportant 
in most central collisions, which result in the overlap of 
the fragment isotopic and isotonic distributions. 

For asymmetric neutron-rich nuclei that have a very 
large N/Z value, more neutrons are pushed to the sur- 
face region and there is a neutron skin structure. In 
the surface region, the neutron density is bigger than the 
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FIG. 4: (Color online) The cross sections of fragments produced in the 140 A MeV 48 Ca + 9 Be and 64 Ni + 9 Be reactions. 
The open circles denote the fragment isotopic a distributions, and the solid ones denote the fragment isotonic a distributions 
of the 48 Ca+ 9 Be reaction. The open triangles denote the fragment isotopic a distributions and the solid triangles denote the 
fragment isotonic a distributions of the 64 Ni + 9 Be reaction. The lines denote the results of the SAA model calculations. 



proton density and the proton-removing will be more dif- 
ficult than the neutron-removing in peripheral reactions. 
The wide spread of the fragment isotopic distributions of 
48 Ca and 64 Ni in the top three panels of Fig. 2] is the ev- 
idence for that. For 48 Ca and 64 Ni, with N/Z equals 1.4 
and 1.286, respectively, there should be a difference in 
their neutron density. The similar shapes of the AA" < 3 
fragment isotonic distributions of 48 Ca and 64 Ni reveal 
the similar proton density distributions in the surface re- 
gions of 48 Ca and 64 Ni, while the shift on the right sides 
of the distributions maybe due to the different diffuseness 
of 48 Ca and 64 Ni according to Ref. [28]. The similarity 
of the left parts of the fragment isotopic distributions re- 
veals the similarity of the neutron density distributions 
in the surface regions of 48 Ca and 64 Ni while the shift 
from the fragment isotopic distributions of 64 Ni to those 
of 48 Ca reveals the difference of the neutron density dis- 
tribution between them. 



IV. SUMMARY 

In summary, by investigating the fragment isotopic and 
isotonic cross section distributions of 140 A MeV 40,48 Ca 
+ 9 Be and 58 > 64 Ni + 9 Be reactions, the evidence of isospin 
effect and its disappearance in projectile fragmentation 
has been found. Some similarity of the fragment isotopic 
and isotonic distributions between not only symmetric 
nuclei but also asymmetric (neutron-rich) nuclei has been 
found. These similarities are related to the similar pro- 
ton and neutron density distributions in the framework 
of the SAA model. These similarities will help to esti- 
mate the fragment production in heavy ion collisions at 
intermediate energy. 
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